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T^t . ABSTRACT 

o: 

Q \ Current data suggest that the central mass densities po and phase-space den- 

"Th \ sities Q = po/o'v ^^ cosmological halos in the present universe are correlated 

^1 with their velocity dispersions ay over a very wide range of ay from less than 

O I 10 to more than 1000 kms~^. Such correlations are an expected consequence 

^ \ of the statistical correlation of the formation epochs of virialized objects in the 

.. \ CDM model with their masses; the smaller-mass halos typically form first and 

. ^ I merge to form larger- mass halos later. We have derived the Q — ay and po ~ (^v 

rS ' correlations for different CDM cosmologies and compared the predicted corre- 

c^ I lations with the observed properties of a sample of low-redshift halos ranging 

in size from dwarf spheroidal galaxies to galaxy clusters. Our predictions are 

generally consistent with the data, with preference for the currently-favored, flat 

ACDM model. Such a comparison serves to test the basic CDM paradigm while 

constraining the background cosmology and the power-spectrum of primordial 

density fluctuations, including larger wavenumbers than have previously been 

constrained. 

Subject headings: cosmology: theory - dark matter - galaxies: formation - galax- 
ies: halos - galaxies: clusters - galaxies: kinematics and dynamics 
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Introduction 



Galaxies and clusters emerge in the CDM model when primordial density fluctuations 
grow to form dark-matter dominated halos in approximate virial equilibrium. For Gaussian- 
random-noise density fluctuations in a cosmologically-expanding, pressure-free gas of col- 
lisionless matter like standard CDM, gravitational instability assembles these halos hier- 
archically. Smaller- mass objects form first, then merge to form larger- mass objects, in a 
continuous sequence of increasing mass. As a result, a statistical correlation between the 
mass and collapse time of the objects which form is expected which depends upon the shape 
and amplitude of the primordial power spectrum P{k) = |5fcp and on the rates of fluctua- 
tion growth in different background cosmologies. This suggests that a measurement of this 
correlation can constrain P{k) and the fundamental cosmological parameters, while testing 
the CDM model. 

The power spectrum is already well constrained at large scales by observations of the 
anisotropy of the Cosmic Microwave Background (CMB) and the mass function of galaxy 
clusters (Bahcall et al. 1999). At large wavenumber k, however, P{k) is still poorly con- 
strained. Recent CMB anisotropy measurements suggest that the universe is flat (e.g. de 
Bernardis et al. 2000, Hanany et al. 2000, Pryke et al. 2001). The measured distance- 
redshift relation for type la SNe (Perlmutter et al. 1999, Riess et al. 2000; Reiss et al. 
2001; Turner and Riess 2001) combines with this to favor a low-density universe with a 
cosmological constant or some other "dark energy" . It is important to establish whether the 
statistical correlation of halo mass and formation epoch predicted at large k by the standard 
CDM model is consistent with these constraints on P{k) at small k and on the fundamental 
cosmological parameters. 

Despite its many successes, the standard CDM model presently faces signiflcant prob- 
lems (e.g. Moore 2001 and refs. therein). For example, CDM simulations flnd singular 
halo density proflles and high abundance of subhalos within halos, in apparent conflict with 
observations, and have difficulty explaining galactic rotation. Recently, modifications of the 
standard CDM model have been proposed to resolve some of these problems, including the 
suggestion of self-interacting dark matter (SIDM) (Spergel & Steinhardt 2000) and other 
CDM variants (e.g. see Wandelt et al. 2000 and refs. therein). In addition to testing the 
standard CDM paradigm while further constraining P{k) at large wavenumber, therefore, 
a measurement of the statistical correlations among the internal properties of halos over a 
wide mass range from dwarf spheroidal (dSph) galaxies (Mhaio ~ 10^ — IO'^Mq) to galaxy 
clusters (Mhaio ~ 10^^ — lO^^M©), expected if halo masses correlate with their formation 
epochs, would sharpen this debate. 

For example, it was recently pointed out (Sellwood 2000, Dalcanton & Hogan 2000, 
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Hogan and Dalcanton 2000, Madsen 2001) that the correlation of the ID velocity dispersion 
(Ty of a halo with the maximum value of its phase space density /max ~ <5 = Po/o-y, where 
Po is the central density of the halo, is a potent diagnostic of the properties of the dark 
matter. As argued in Sellwood (2000), Liouville's theorem applied to collisionless dark matter 
guarantees that processes like baryonic dissipation will not change /max significantly. In that 
case, /max can be used as a diagnostic even for halos that are no longer dark-matter dominated 
in the central region following baryon cooling and compression. Of related interest, there 
are recent claims that the core densities of all galaxy and cluster halos are approximately 
the same, independent of the mass of the halo (Firmani et al. 2000, Avila-Reese et al. 2000, 
Kaplinghat, Knox, & Turner 2000). A central density which is the same for halos of different 
mass would impose strong constraints on the nature of the dark matter, if this hypothesis is 
correct. 

In what follows, we predict correlations of the maximum phase-space densities Q and of 
the central densities po of cosmological halos with their velocity dispersions ay, for different 
CDM models. We will show that these predicted correlations are consistent with the observed 
properties of halos spanning a very large range of mass scales from dSph galaxies to galaxy 
clusters, and that the current data has the power to constrain the background universe and 
power spectrum of the CDM model. Finally, we shall discuss the hypothesis that halos have 
a core density that is independent of their mass. 



2. The Q — ay and po — ay Correlations for Cosmological Halos Observed Today 

We have developed an analytical model for the postcollapse equilibrium structure of 
virialized objects which condense out of a cosmological background universe, either matter- 
dominated or flat with a cosmological constant (Shapiro, Iliev & Raga 1999, "Paper I" ; Iliev & 
Shapiro 2001b, "Paper II"). The model is based upon the assumption that cosmological halos 
form from the collapse and virialization of "top-hat" density perturbations and are spherical, 
isotropic, and isothermal. This leads to a unique, nonsingular, truncated isothermal sphere 
(TIS), a particular solution of the Lane-Emden equation (suitably modified when A 7^ 0). 
The size rt and velocity dispersion ay are unique functions of the mass and redshift of 
formation of the object for a given background universe. Our TIS density profile flattens 
to a constant central value, Pq, which is roughly proportional to the critical density of the 
universe at the epoch of collapse, with a small core radius Tq ~ Ti/30 (where ay = AnGpoVQ 
and To = rKing/3, for the "King radius" rxing, defined by Binney and Tremaine, 1987, p. 
228). 

This TIS model reproduces many of the average properties of the halos in CDM N- 
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body simulations quite well, suggesting that it is a useful approximation for the halos which 
result from more realistic initial conditions. We have elsewhere compared this model with 
simulation results and applied it successfully in a variety of ways, including the derivation 
of the mass-radius-temperature virial relations and integrated mass profiles of X-ray clus- 
ters deduced from numerical CDM simulations by Evrard, Metzler & Navarro (1996) and 
Mathiesen & Evrard (2001) (Papers I, II), the rotation curves of dark-matter dominated 
dwarf galaxies and the observed correlation between their maximum velocity, v^ax, and the 
radius, r^-a^x, at which it occurs (Iliev & Shapiro 2001a), and the mass model of cluster 
CL 0024 determined by strong gravitational lensing measurements (Shapiro & Iliev 2000). 
The TIS mass profile agrees well with the fit to N-body simulations by Navarro, Frenk & 
White (1996; "NFW") (i.e. fractional deviation of ~ 20% or less) at all radii outside of a 
few TIS core radii (i.e. outside a King radius or so), for NFW concentration parameters 
4 < cnfw ^7. As a result, the TIS central density po can be used to characterize the core 
density of cosmological halos, even if the latter have singular profiles like that of NFW, as 
long as we interpret po, in that case, as an average over the innermost region. 

The TIS model yields the central density po, velocity dispersion ay, size rt and core 
radius tq for any given halo as unique functions of its mass M and collapse redshift Zcow- This 
defines a "cosmic virial plane" in the three-dimensional space of halo parameters {po,rQ,(Tv) 
[i.e. 2lnav — Inpo ~ 21nro — ln(47rG') = 0], in the terminology of Burstein et al. (1997), 
and determines the halo size, mass and collapse redshift for each point on this infinite plane, 
the ensemble of all possible virialized objects. In hierarchical models of structure formation 
like CDM, however, the statistical correlation of halo mass and collapse redshift determines 
the expected distribution of points on this cosmic virial plane. We shall combine our TIS 
model with the well-known Press-Schechter (PS) approximation (Press & Schechter 1974) 
for Zco\\{M) - the typical collapse epoch for a halo of mass M - to derive the statistical 
distribution for the Q — cry and po ~ ^v correlations, as follows. 

According to PS, the fraction of matter in the universe which is condensed into objects 
of mass > M at a given epoch is erfc (z//2-^/^), where u = 5cvit/o'{M), a{M) is the standard 
deviation of the density fluctuations at that epoch, according to linear theory, when filtered 
on mass scale M, and (^crit is the fractional overdensity of a top-hat perturbation when this 
linear theory is extrapolated to the time of infinite collapse in the exact nonlinear solution. 
The "typical" collapse epoch for a given mass is that for which a{M) = 6crit (i-e. z/ = 1). 
For a given 2;coii, this defines a typical mass scale: M^ = M{v = 1). 

For any z/, the PS prescription above yields a unique Zcou for halos of a given mass, thus 
completely determining the Q — ay and po — cry correlations. The TIS-I-PS predictions for the 
Q — ay correlation for halos observed today are shown in Figure 1, for u = 0.5 — 3. For any 
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z/, there is some value of ay = avfl for which Zcow = 0. For fluctuations which become halos 
observed today with ay > ay^, the most hkely collapse epoch is assumed to be Zcow = 0. 
Hence, curves of constant u merge with the straight line for ZcoW = for this range of ay 
values. 

For comparison, the observed properties of a sample of low-redshift, dark- matter domi- 
nated halos of galaxies and clusters are also plotted in Figure 1, using data from the following 
sources. Kormendy and Freeman (1996; 2001) have compiled and interpreted mass models 
fitted to galaxy rotation curves for 49 late-type spirals of type Sc-Im and velocity dispersions 
and core radii for 7 dSph galaxies, chosen to minimize the need to correct the deduced halo 
parameters for the effects of baryonic dissipation. Data for one additional dSph galaxy, Leo 
I, is from Mateo et al. (1998). For low-redshift clusters, we have adopted central total mass 
densities which are based upon the central gas densities and baryonic-gas-mass fractions de- 
rived from X-ray brightness profile fits for 28 nearby clusters by Mohr, Mathiesen & Evrard 
(1999), for which velocity dispersions were tabulated by Girardi et al. (1998) (26 clusters) 
and Jones & Forman (1999) (2 clusters). The galaxy and cluster halo profiles adopted by 
these authors in fitting data to derive their central densities are very similar to the TIS halo 
profile we have used to make our predictions, so these central densities are suitable for direct 
comparison. 

The agreement in Figure 1 between our TIS+PS model predictions and the observed 
halos is excellent for the entire range of halo masses from dSph galaxies to galaxy clusters. 
The observed scatter of the data is natural from a theoretical point of view, due to the 
Gaussian random statistics of the initial density fiuctuations. A comparison of results for 
different models suggests that the currently-favored, flat ACDM model {Qq = 0.3, Aq = 
0.7, h = 0.65) is preferred. The COBE-normalized, low-density, matter-dominated OCDM 
model {Qq = 0.3, Xq = 0,h = 0.65), with a primordial power spectrum index of Up = 1.3, a 
tilt which allows it to match the galaxy cluster abundance aX z = (e.g. Wang et al. 2000), 
requires that most galaxy halos observed today be low-z/ (i.e. late-collapsing), rather than 
typical. This may indicate that P{k) for this model is too large on galactic mass scales. The 
untilted OCDM model (up = 1), which predicts an incorrect cluster abundance ai z = 
(i.e. not "cluster- normalized" ) , requires, instead, that all halos observed today be high-z/, 
corresponding to rare, high-density peaks of the primordial density fleld. In that case, P{k) 
may be too small on cluster mass scales and below. Finally, the cluster-normalized Einstein- 
de Sitter (EdS) model ("SCDM") shows good agreement with this halo data on galaxy scales, 
but there are many cluster data points for which the Q-values fall below the limiting values 
predicted to correspond to z^oii = 0. 

In Figure 2, we compare the TIS+PS prediction for the po — <Jy correlation with the same 
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halo data. Once again, the data are generally consistent with the theoretical prediction, with 
the same trends for different CDM models as identified above for the Q — ay correlation. 
Hence, the flat, ACDM model is preferred. The data are clearly inconsistent, however, with 
the constant-core-density hypothesis, since observed core-densities show a clear downward 
trend from dSph galaxies to larger galaxies and clusters, decreasing by more than three 
orders of magnitude as ay increases by two orders of magnitude. The appearance of a core- 
density that is approximately independent of ay at large values of ay (i.e. from massive 
galaxy to cluster scales) is easily understood in the TIS+PS model if the Gaussian statistics 
of fluctuation amplitudes are considered. Since large galaxies and galaxy clusters observed 
today formed from high-z/ fluctuations, the most probable time for them to collapse is close 
to the present. Therefore, our model predicts them all to have similar central densities, since 
Po is proportional to the mean density of the universe at the collapse epoch. 

This success of the TIS+PS model in predicting the observed correlations can be under- 
stood by an analytical argument, as follows. According to Paper II, the exact dependences of 
the TIS model core densities and velocity dispersions on M and Zcow are well approximated 
by Po = 1.80 X 10^[F(2eoii]'Vcrit,o and 4 = 1.1 x lQ\M/lQ^^MQh-^f/^F{z,,x{)]mv^s-\ 
where Pcrit.o = SH^/{87rG) and the function F{zcoii) is deflned by equation (84) of Paper 
II. The quantity F^ is just the density of a top-hat perturbation after it collapses and viri- 
alizes in the standard uniform sphere (SUS) approximation (as measured in units of the 
EdS value for Zcou = 0, ISvr^pcrit.o)- For the EdS case, F = (1 + Zcow), while for the low- 
density, matter-dominated and flat cases, F -^ Qq (1 + Zcon) for early collapse. If we 
approximate the power-spectrum of density fluctuations at high redshift (e.g. just after re- 
combination) as a power-law in wavenumber k, P{k) oc A;", and deflne a mass M oc k~^, 
then F oc M-(3+")/^ if n = n^s = -3(2^^ + 1), where yp = (rflnF/dlnM)exact at the rele- 
vant mass scale. For all masses in the EdS case and for masses which collapse early in the 
low-density, matter-dominated and flat, universes, yp reduces to yo- = ((ilncr/(ilnM)exact, 
and (1 + -Zcou) oc a{M) oc M"*^^"*""-'/®, where a{M) is evaluated at the same cosmic time for 
all masses. The dependence of Ues and Zcow on M for 1-a fluctuations is shown in Figure 3 
for ACDM, along with the approximate Ucs which results if yp is replaced by ycr, which 
shows that the latter is a very good approximation for all masses M < 10^"^ MqH^^ . In terms 
of this power-law model, the TIS model then yields po oc M"*^"^^^/^, ay oc M*^^"")/-*^^, and 
Q oc M~^"'~^'^''*, which combine to give 

Q = Q,{ay/ay,,r-\ (1) 

and 

Po = po,*(o-y/o-y,^)", (2) 

where a = 6(n + 3)/(n — 1), except that a = for M > M^{z = 0), for which z^oii = is 
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assumed. 

For the currently-favored, flat, ACDM model, for example, Ucs ~ —2.6, —2.4, —2.2 for 
M = M^ = 10\W^,W^MQh-\ respectively, and, hence, a ^ -0.7,-1.1,-1.8 for these 
masses, while ay,* = 4.5, 33, 114 km s-\ Q^ = 1.4x10^3^ 5.7x 10"^ 2.4x 10"^ M0pc-3(km s'^)-^, 
and po,* = 0.12, 0.02, 0.0035 M0pc-^ respectively. For all masses M > M^{z = 0) ^ 
IO^^H^^Mq, our model predicts po(-2coii = 0) = 0.00143 MqPc~^, and (^(-Zcoii = 0) = po(-2coii = 

o)K. 

Dalcanton and Hogan (2000) predict that Q oc ay^ for all ay, which in our model 
is achieved only for M — > 0, for which rics = —3, and in the high-mass limit, for which 
-Zcoii ~ typically. Similarly, according to equation (2), the TIS + PS model predicts that 
the core-density is the same for halos with different ay only in the small-mass limit in which 
Ues = —3 and the high-mass limit, for which ZcoU ~ 0, with a large difference between the 
values of po in these two limits. The generally good agreement between this prediction and 
the data, as plotted in Figure 2, argues against the hypothesis that all halos share a single, 
universal core density. The putative universal core density may have escaped detection by 
the observations quoted here if it is confined to a region interior to the smallest radii probed 
by those observations, but it must be at least as large as the highest central densities reported 
in Figure 2. These high central densities are well in excess, of the universal core density of 
O.O2M0PC-3 claimed by Firmani et al. (2000). 

In the future, as more and better data become available, it should be possible to refine 
the comparison between the correlations predicted here and the observed properties of galax- 
ies and clusters. We have demonstrated that such a correlation is expected, that current 
data are consistent with the predicted correlations for the CDM model, with preference for 
the currently-favored ACDM model over CDM in other background cosmologies, and that 
no significant deviation of the primordial power-spectrum shape from the scale-invariant, 
Harrison-Zeldovich shape {up = 1), the standard prediction of inflationary cosmology, is 
required. In addition, we have demonstrated that the data do not support the idea that all 
halos have the same core density, independent of their mass or formation epoch. 

We thank Hugo Martel, Sandra Faber and Ken Freeman for discussion and John Kor- 
mendy for providing his tabulation of data points for galactic halos used here. This research 
was supported by grants NASA ATP NAG5-7363, -7821, and -10825, TARP 3658-0624- 
1999, European Community RTN contract HPRN-CT2000-00126 RG29185, and NSF INT- 
0003682. 
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Fig. 1. — Maximum phase-space density Q = po/uy versus velocity dispersion ay for halos 
observed today, as predicted for various CDM universes by the TIS + PS model (solid 
curves) for halos formed from v-a fluctuations, as labelled with the values of z/, for v = 
0.5, 1, 2, and 3. In each panel, lines representing halos of different mass which collapse 
at the same redshift are shown for the case Zcow = 0, as labelled. Each panel represents 
different assumptions for the background universe and primordial density fluctuations, as 
labelled: COBE-normalized ACDM (Aq = 0.7, Qq = 0.3) (upper left), cluster-normalized 
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SCDM (fio = 1) (upper right), and COBE-normalized OCDM {Qq = 0.3) (lower panels), all 
assuming h = 0.65 and primordial power spectrum index n^ =1 (i.e. untilted), except for 
OCDMi.3, for which Up =1.3. Data points represent observed galaxies and clusters, taken 
from the following sources: (1) galaxies from Kormendy & Freeman (1996; 2001) (open 
triangles); (2) Local Group dwarf galaxy Leo I from Mateo et al. (1998) (filled square); (3) 
cluster velocity dispersions from Girardi et al. (1998) and Jones & Forman (1999), central 
densities from Mohr et al. (1999) (crosses). 

Fig. 2. — Same as Figure 1, but for central density po versus velocity dispersion cy. 

Fig. 3. — (a) (upper) Effective logarithmic slope, Ues, of the density fiuctuation power spec- 
trum versus halo mass M = M^ (i.e. for 1-a fluctuations) for ACDM, based upon yp (solid 
curve) and the approximation which uses y„, instead (dashed), (b) (lower) Typical collapse 
redshift Zcow for halos of mass M* in ACDM. 
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